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STUDY  OF  THE  EFFECT  OF  LOCAL  GEOPHTSICAL  PHEHOMEMA  OI^IfmEWCES 
LATITUPy-S  awp  MTAV  LATITUDES  OF  STATIOHg  LOCATED  IH  Tffi; 
VTCINITY  of  a  common  MERIDIAH  in  the  period  1 962*1  fr?5 

Barbara  Kolaczek,  Ronan  Galas,  Marcin  Barllk,  Magdalana  Dukwics 

I.  Introduction 

Changes  of  geographical  latitude  are  caused  mainly  by 
free  or  forced  nutation  of  the  axle  of  Earth's  rotation.  In 
addition,  changes  of  geographical  latitude  are  affected  by 
a  number  of  other  factors  such  as  errors  in  astronomical  constants 
and  in  declination  of  stars,  instrumental  errors,  Inaccurate 
theory  of  nutation  and  refraction,  and  also  local  movements  of 
Earth's  cruet  and  changes  in  the  plumb  direction. 

Local  movements  of  Earth's  crust  and  changes  in  the  plumb 
(vertical)  direction  result  mainly  in  secular  and  irregular 
changes  Of  the  mean  geographical  latitude,  tliat  is  the  latitude 
freed  from  periodical  variations.  The  mean  latitude  of  a  station 
is  also  influenced  by  irregular  and  secular  instrumental  changes 
and  refractional  anomalies.  On  one  hand,  local  changes  of 
latitude  are  superimposed  on  changes  caused  by  movements  of 
poles,  complicating  the  picture.  On  the  other  hand,  local  changes 
of  latitude  contain  information  on  the  local  movements  of  the 
Earth's  crust  -  both  those  connected  with  tectonic  plate  displace- 

1 

ments  and  wit3i  local  mass  movements  caused,  for  instance,  by 
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earthquakes.  Investigation  of  the  local  changes  of  latitude  Bay, 
therefore,  facilitate  studies  of  the  BOveBents  of  poles  and  Bay 
contrlbv^e  to  studies  of  tectonic  plate  displaceaents.  The  ualfoni 
observations  of  changes  In  latitude  over  a  period  of  years 
represent  a  unique  bank  of  data  for  goodynaalc  studies. 

The  secular  local  changes  of  latitude  for  sobs  stations 
have  been  determined  (/V,  /5/»  /6/,  /H/t  /^7/,  /21/)  and  they 
are  usually  Interpreted  as  Boveoents  of  tectonic  plates.  Irregular 
sudden  changes  of  the  mean  latitude  of  the  order  of  091  (several 
months)  have  also  been  noted  A/,  /1 3/*  However,  the  reasons  for 
their  appearance  remain  unknown.  Investigation  of  correlation 
between  the  local  changes  of  latitude  and  the  better-known  now 
extents  of  displacement  of  tectonic  plates,  the  occurrence  of 
earthquakes  near  stations  located  in  active  seismic  regions, 
and  finally  changes  in  plumb  line,  would  enable  to  explain  Bain 
reasons  for  local  changes  of  latitude. 

To  explain  the  character  of  the  local  changes  of  latitude 
we  analysed  in  this  work  changes  of  mean  latitude  and  differences 
between  mean  latitudes  for  stations  located  near  the  Beridlan 
of  station  Jozefoslaw.  The  placement  of  these  stations  is 
particularly  suitable  for  this  investigation.  Moreover,  this 
work  covers  also  experimental  results  of  graviBetric  measure¬ 
ments,  carried  out  at  the  meridional  base  Jozefoslaw  to  deterBine 
meridional  component  changes  in  the  plumb  line  and  their  effect 
on  changes  of  latitude. 


TT  r!hany«t«  of  Mean  latltndo  and  dlffT»ncee  of  »»aa  latltudae 
^n-r  Btatlone  located  near  »  common  MTldlan 

Differences  of  geographical  latitude  and  nean  latitude 
for  stations  located  In  the  vicinity  of  a  conmon  meridian  are 
particularly  suitable  for  studies  of  the  effect  of  local  factors 
since  the  effect  of  the  movements  of  pole  Is  here  eliminated. 
Changes  In  differences  of  the  mean  latitude  of  Jozefoslam  and 
Belgrade  have  been  noted  for  the  first  time  by  Telekl  /20/, 

In  the  present  work  we  analyzed  changes  of  mean  latitudes  and 
their  differences  for  four  stations  located  near  a  common 
meridian,  namely,  Turku,  Jozefoslaw,  Belgrade  and  Pecny, 

Turku  is  situated  in  the  area  of  large  vertical  movements 
of  the  Northern  shores  of  Baltic,  while  Belgrade  is  located 
In  a  eeismlcally  active  region,  near  the  boundaries  of  two 
tectonic  plates.  In  these  areas,  local  movements  of  Earth's 
crust,  and  the  associated  changes  in  latitudes  of  stations, 
can  be  expected.  For  cov.parlson,  changes  of  mean  latitudes  of 
several  European  stations  have  also  been  calculated.  Table  1 
gives  .coordinates  of  the  considered  stations. 

On  the  basis  of  data  published  by  IPMS  in  "Annual  Report" 
for  the  period  1962-1975  /9/  and  in  "Monthly  Notes"  for  the 
period  1973-1975  /10/  we  calculated  mean  latitudes  of  these 
stations,  using  an  0rlov-0R2  filter.  By  application  of  the 
method  of  harmonic  analysis  we  calculated  also  the  so-called 
filtered  latitude  of  these  stations  /8/,  The  filtered  latitude 


was  obtained  ae  a  reault  of  ellninatlon  froa  the  observational 
data  of  all  the  latitudes  In  periodic  vslyhs  whose  amplitude 
exceeded  the  mean  error  of  determination,  that  is  the  value  of 
0?004-0l*008.  The  determined  periodic  vmlmhm  of  the  variations 
of  latitude  for  the  considered  stations  in  the  period  1962-1975 
are  presented  in  Table  2,  Graphs  of  the  mean  latitudes  counted 
with  the  Orlov  filter  and  of  the  filtered  latitudes  are  shown 
in  Figure  1 .  On  Figure  1  we  marked  also  by  arrows  the  moments 
of  occurrence  of  larger  (M  >5.5)  earthquakes  in  the  region  of 
Belgrade  (at  a  distance  shorter  than  2000  km),  and  the  instrumental 
changes  reported  by  observatories  in  Turku  and  Belgrade. 

we  calculated  also  differences  of  latitudes  as  well  as 

differences  of  mean  latitudes  calculated  with  Orlov  formula 
for  the  considered  stations.  The  graphs  of  these  latter 
differences  are  shwon  in  Fig.  2.  Changes  of  differences  in  the 

mean  latitudes  represent  now  only  the  relative  changes  of  the 

mean  latitudes  of  stations  caused  by  the  local  factors. 

Variations  of  mean  latitudes,  counted  with  Orlov  filter, 
and  their  differences  for  the  remaining  European  stations  are 
shown  in  Figures  3  and  4. 

The  analysis  of  variations  in  mean  latitudes  and  their 
differences  enables  to  draw  the  following  conclusions. 

1  •  The  largest  variations  of  the  mean  latitude  occur  in 
Belgrade,  that  is  the  station  situated  in  a  seismlcally  active 
zone.  On  Figure  1b  one  can  notice  a  correlation  between  the 
sudden  decrease  of  the  mean  latitude  of  Belgrade  and  the 


earthquake  in  Skopie  No.  0  (Table  3)f  as  well  as  other  saaller 
earthquakes  No.  2,  4t  7he  earthquake  No.  12,  whose  epicenter 
was  located  closer  to  Belgrade  than  all  the  renalning  earth" 
quakes,  coincides  also  with  a  sudden  change  of  the  mean  latitudle 
and  also  with  the  then  introduced  thermal  protection  of  "libel*. 
The  seismic  activity  in  the  region  of  Belgrade  meridian  was 
more  pronounces  in  the  years  1964-1966,  1968-1970  and  1972 
and  it  correlated  with  larger  changes  of  the  mean  latitude  of 
Belgrade,  as  well  as  differences  of  mean  latitude  of  Belgrade 
and  Jozefoslaw. 

2.  Variations  of  the  mean  latitude  of  stations  as  well  as 
their  differences  are  smaller  in  the  case  of  continental 
stations  (Jozefoslaw,  Pecny,  Paris),  and  larger  for  stations 
located  nearer  to  seashores  (Turku,  Hamburg,  Hccle). 

3.  On  the  basis  of  available  data  it  is  difficult  to  explain 
the  reasons  for  such  sudden  changes  in  latitude  as  for  Turku 
In  1972,  for  Vccle  in  1965  and  1972,  and  for  Hamburg  in  1972. 

It  is  known  that  Turku  lies  in  the  region  of  large  vertical 
movements  of  seashore  of  the  Baltic  Sea,  hence  in  an  unstable 
area  where  local  displacements  of  mass  in  the  region  of  station 
can  take  place. 

Because  of  the  short  period  of  time  and  special  choice  of 
stations,  the  analyzed  variations  of  the  mean  latitude  did  not 
allow  to  draw  conclusions  with  respect  to  the  displacements 
of  tectonic  plates. 


4.  The  curves  of  variations  of  the  filtered  latitude  of 
station  (Figure  1 )  enable  to  register  in  aore  detail  the  irre¬ 
gular  changes  than  does  the  mean  latitude  counted  with  very 
approximating  (smoothing)  Orlov  filter,  as  can  be  seen,  for 
instance,  in  the  case  of  Turku.  Moreover,  spectrally  filtered 
curve  of  changes  in  latitude  provides  an  information  about 
secular  and  irregular  changes  of  latitude  throughout  the  whole 
period  considered,  without  omitting  the  one-year-and-a-half 

or  three-year  end  periods,  as  happens  in  the  cases  of  the  Orlov 
or  BIH  filters.  This  fact  is  of  prime  importance  in  studies  of 
variations  of  the  mean  latitude. 

5,  In  addition  to  main  parameters  of  periodic  expressions 
for  changes  of  latitude,  that  is  the  Chandler  expression  and 
the  annual  and  semiannual  expressions,  the  parameters  of  many 
other  periodic  expressions  were  determined  as  a  result  of 
harmonic  analysis  (Table  2),  Some  of  them  are  known  in  the 
literature  and  some  attempts  of  physical  interpretation  have 
been  made  for  them,  as  noted  in  Remarks  of  Table  2.  However, 
in  the  interpretation  of  these  expressions  one  has  to  remember 
that  some  artificial  periodic  expressions  can  arise  as  a  conse¬ 
quence  of  methods  of  treatment  of  experimental  data  and  of  the 
application  of  the  method  of  harmonic  analysis  to  a  process 
which  is  not  exactly  fully  suitable  for  it.  In  values  of  the 
mean  latitude  counted  with  the  aid  of  Orlov  formula  for  stations 
Jozefoslaw,  Turku,  Belgrade  and  Pecny  there  occurred  quite 


dletinct  long-period  expressions  having  the  periods  2, 6-2, 8 
years  and  the  amplitudes  0y01-0'.'02  (/1 8/  »  /2/), 

Further  analysis  of  variations  of  the  mean  latitude  of 
satatlons  considered  here  and  of  other  stations  located  in 
different  geophysical  regions  may  help  to  explain  the  character 
of  local  changes  in  the  latitudes  of  stations,  and  also  may 
eventually  give  selection  of  latitudinal  stations  whose 
observational  data  would  not  contain  large  local  disturbances, 

HI,  Determination  of  c  ^nxes  in  the  plumb  line  at  the  meridian 
of  Jozefoslaw  through  gravimetric  measurements 

Ill.t,  Introduction 

Lunisolar  influences  and  dislocation  of  mass  inside  the 
Earth,  including  the  displacement  of  tectonic  plates,  exert 
effect  upon  the  orientation  of  equipotentlal  surface  at  the 
point  of  observation,  that  is  on  the  orientation  of  the  actual 
plumb  line.  Changes  in  the  direction  of  pliimb  line  di  enter 
directly  into  determination  of  the  astronomical  geographical 
latitude,  since  ^  ,  To  present  mathematically  the  form 

of  this  dependence  we  shall  make  use  of  the  well  known  equation 
describing  the  change  of  the  plumb  line  on  account  of  lunar 
and  solar>  effects 


— ~Ak - —  un2z, 

2  g  R* 


(1) 
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where  -  gravitational  constant,  yn  -  mass  of  Sun  (Moon), 

Y  -  radius  of  Earth,  -  distance  to  Sun  (Moon),  Z  -  zenith 
distance  of  Sun  (Moon),  where  %  and  L  denote  Love's 

numbers  characterizing  the  compressibility  of  Earth,  The  value 
of  A  'I.13*  and  in  our  considerations  we  shall  assume  it  to  be 
a  constant  with  time. 

Instead  of  the  potential  of  attraction  of  Sun  (Moon),  one 
can  use  in  the  above  equation  the  potential  of  attraction  of  an 
additional  mass  inside  the  Earth.  It  could  le  located  under  the 
surface  periodically  as,  for  instance,  ground  water,  or  it  could 
anpear  in  a  sudden  manner  as  it  happens  in  the  case  of  an 
earthquake.  This  movement  should  reflect  itself  also  in  changes 
of  equlpotential  surface  and  direction  of  the  plumb  line,  provided 
these  displacements  involve  sufficiently  large  masses.  The  directtn 
of  the  plumb  line  can  be  calculated  from  the  equations  of 
Venlng  Meinesz,  The  meridional  component  4^  of  the  deviation  of 
plumb  line  with  gradient  of  gravimetric  anomaly  /^^can  be 
represented  by  equation 


di  -  0"105  d-  =  V:(a25{Ag,-AgM), 


(2) 


where  ^  denotes  the  distance  of  points  /V  and  S  of  the  base, 
and  denotes  Faye's  gravimetric  anomaly  determined  at  those 
points.  The  coefficient  0,105  arises  from  theory  of  the  method 
of  Vening  Meinesz, 
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After  time  Interval  between  consecutive  observations 
and  ,  changes  in  gravimetric  anomaly  will  cause  changes 
in  the  direction  of  the  plumb  line  by  the  amount 


Ji  -  a'os25(4r.-4r»). 


The  above  considerations  are  valid  only  for  the  model  of 
rigid  Earth.  For  the  real  Earth  the  last  equation  should  be 


modified  by  introduction  of  the  Love's  numbers.  The  measured 
differences  of  acceleration  are  Crs  *  1  .20  times 

larger  than  those  for  the  model  of  rt.gid  Earth,  In  turn,  the 
observed  values  of  geographical  latitude  are 
larger  than  those  obtained  from  the  equation  (3).  Finally,  we 
get  the  working  formula  for  the  described  method: 

Jf-J{»Or0523-^  Wf  »0:'0f94W,.  (4) 


A  simple  analysis  of  this  equation  will  show  that, in  order 
to  detect  a  change  of  the  mean  latitude  of  the  observatory 
of  the  order  of  0V001 ,  the  error  of  determination  of  changes 
in  acceleration  should  not  exceed  20  p  gal.  The  required 
accuracy  of  the  determination  of  dictates  and  limits 

the  method  of  gravimetric  work  and  the  type  of  used  instruments. 
This  accuracy  can  be  achieved  by  means  of  modern,  astatlslzed 
quart?  instruments. 
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III. 2.  Results  of  determination  of  changes  in  the  direction  of 
plumb  line  at  Jozefoslaw  in  the  period  of  one  year 

On  the  basis  of  given  above  considerations  ve  performed 
experimental  determination  of  changes  of  the  mean  geographical 
latitude  of  Jozefoslaw  caused  by  movements  of  the  direction  of 
plumb  line  (vertical).  For  this  purpose,  we  established 
a  meridional  gravimetric  base  containing  seven  measurement 
points,  namely: 

NI4  -  a  point  14  km  North  of  Jozefoslaw  station;  N6  - 
a  point  6  km  North  of  station;  N3  -  a  point  distant  3  km  North 
of  the  station;  and  similarly  points  SI 4  >  S6  and  S3  located 
correspondingly  at  the  distance  of  14»  6  and  3  kilometers  South 
of  the  Jozefoslaw  station. 

All  these  points  were  chosen  following  the  strictest  require¬ 
ments  which  apply  to  establishing  a  network  of  basic  gravimetric 
points.  The  distances  between  points  of  the  base  were  such  that 
by  means  of  repeated  gravimetric  measurements  one  could  detect 
dislocations  of  mass  near  the  Jozefoslaw  station  -  point  0, 
movements  of  mass  of  the  Earth  crust,  and  movements  of  even 
deeper  layers  of  Earth  to  Mohorovici^  boundary  (points  NI4  and 
si  4).  Feitelson  et  al,  /7/  suggested  the  utilization  of 
a  surface  gravimetric  picture  of  area  up  to  200  km  from  the 
station  for  this  purpose.  Such  a  large  outline  of  the  picture 
does  not  allow  to  carry  out  fast  periodic  repeats  of  measurements 
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to  determine  changes  in  the  direction  of  plumb  line.  Therefore, 
for  this  work  we  adopted  the  concept  of  measurements  on 
neridional  gravimetric  base,  so  much  more  that  we  are  primarily 
interested  in  changes  of  the  meridional  component  of  the 
direction  of  plumb  line.  Measurements  of  differences  of  the 
acceleration  at  the  base  were  performed  by  means  of  gravimeters 
Worden  No,  867,  or  Sharpe  Nos,  129,  185,  202  and  203* 

The  uniform  scale  of  determined  values  was  ensured 
through  standardizing  gravimetry  by  the  method  of  inclination 
on  eraminator.  The  measurements  carried  out  so  far  were  in  four 
series:  10-11  June  1976,  6-9  November  1976,  9  March  1977,  and 

5.7  Mav  1977,  In  each  series  the  observations  were  made  in  three 
rounds,  and  the  spread  of  measurements  over  several  days  was 
intended  to  provide  a  better  averaging  of  the  results  of 
observations.  Appropriate  instrumental  adjustments  and  lunlsolar 
corrections  were  introduced  into  the  results.  Table  4  presents 
differences  AJ  ,  and  their  mean  errors  for  particular  sections 
of  the  base  and  for  particular  series  of  measurements.  The  Table 
lists  also  changes  in  the  direction  of  plumb  line  (vertical) 
calculated  according  to  the  formula  (4). 

AB  a  result  we  obtained  the  following  changes; 

In  the  period  from  June  1976  to  November  1976 

=  -  0,”0024  1  0','0015; 

In  the  period  from  November  1976  to  March  1977 


=  +  oyooi9  +  oyooi8; 


In  the  period  fro®  March  to  May  I977 

•  0V0030  +  050012, 

The  obtained  reeults  Indicate  the  possibility  of  registra¬ 
tion  by  this  method  of  changes  In  geographical  latitude  of  the 
order  of  0501  caused  by  movements  of  the  plumb  line.  These 
first  experimental  results  do  not  Justify  yet  physical 
interpretation  of  the  registered  changes  in  mean  latitude. 
Small  changes  of  the  mean  latitude  require  the  introduction  of 
a  large  number  of  measurements  in  short  Intervals  of  time.  The 
uee  of  a  larger  number  of  instruments  mould  allow  to  eliminate 
systematic  errors  of  the  instruments.  Performing  of  such 
gravimetric  measurements  at  latithdlaalj^tatlons  situated  in 
seismically  active  regions,  as  for  instance  Belgrade,  would 
allow  a  deeper  geophysical  interpretation  of  changes  in  the 
mean  geographical  latitude. 

Table  1.  Coordinates  of  the  considered  latitudinal 
stations . 


T 

X 

9 

Turku 

wrr 

Uccle 

-0*17" 

50*48' 

lAacfoshkW 

52*06' 

Paris 

-0*09" 

48*50’ 

Belgmd 

-1*22" 

44*48' 

Neuchatel 

-0*28" 

47*04' 

PeCny 

-irsB" 

4y34' 

Hamburg 

-0*40" 

53*36' 

Key:  (1)  Station. 
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Table  2.  Periodicities  of  changes  in  the  mean  latitude  of 
stations  located  at  a  meridian  common  to  station  Jozefoslaw» 
determined  in  the  process  of  filtration. 


Jtefodiw 

P  A 

‘hlffctt 

P  A 

Pefay 

P  A 

Belgrad 

P  A 

UwaciO) 

1:193 

ori3i 

inw 

orii2 

i;i80 

orio9 

i:i88 

0:'143 

w.  Chanlenta.) 

1.00s 

0.162 

1.0m 

0.101 

0.996 

0.125 

1.030 

0.124 

w.  rocznyO) 

0.609 

0.029 

0.5012 

0.031 

0.595 

0.010 

0.530 

0.034 

w. 

— 

_ 

0.330 

0.047 

0.321 

0.034  1 

— 

— 

— 

— 

.330 

0.017 

.338 

.014  1 

1  wplyw  anoimlii 

— 

— 

— 

— 

.394 

.023 

.375 

.018  I 

1  nfrakcji  (37) 

0.427 

0.013 

0.406 

0.011 

.424 

.015 

.415 

0.25 

— 

— 

.443 

.010 

.439 

.021 

.424 

0.14 

.473 

.013 

— 

— 

.463 

.014 

.436 

.028 

— 

— 

— 

.487 

.027 

.468 

.019 

wplyw  okolodo- 

— 

.313 

.011 

— 

— 

.504 

.021 

bowej  nutacji 

.5«8 

.011 

.385 

.008 

.551 

.021 

.536 

.015 

swobodnej  [11]. 

— 

— 

— 

— 

— 

— 

.577 

.014 

[I2J.  fl3J  (^) 

— 

— 

— 

— 

.632 

.013 

.624 

.014 

.092 

.015 

.680 

.013 

.638 

.019 

.638 

.012 

.767 

.016 

.723 

.010 

.709 

.011 

.784 

.016 

.909 

.010 

.882 

0.16 

.883 

.022 

.881 

0.17 

1 

.930 

.031 

— 

— 

__ 

.962 

.055  1 

wplyw  bt^du  sta-  1 

1.0S0 

.016 

- 

1.129 

.017 

1.075 

.017  i 

lej  nutacji  |18)^J  ■ 

1.312 

.027 

1.324 

.021 

1.334 

.029 

1.268 

.034  j 

wplyw  oscylacji  | 
oceanu  [16]  ( 

1.316 

.016 

— 

1.448 

.01.1 

1 

-- 

— 

— 

1.650 

.013 

1.637 

.011 

[18) 

1.9.19 

.014 

1.813 

.012 

2.020 

.017 

2.090 

.011  1 

wplywy  nictcorolo- 1 

2.444 

.009 

2.374 

,034 

2.507 

.015 

— 

-  ( 

gicinc 

Key:  (1)  Remarks;  (2)  Chanler  period; 

(3)  Annual  period;  (4)  Semiannual 

period;  (5)  Effect  of  refraction  anomaly  [37];  (^)  Effect 
of  diurnal  free  nutation  [11],  [12],  [15];  (7)  Effect  of 
the  error  of  constant  nutation  [I8];  (8)  Effect  of 
oscillation  of  ocean  [I6];  (9)  Meteorological  effects  [2]. 
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Figure  lb.  Curves  of  changes  in  the  mean  and 
filtered  latitude  of  stations  Peony  and  Belgrade. 
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Figure  2.  Curves  of  changes  in  the  mean 
latitude  of  stations  located  in  the  vicinity 
of  a  common  meridian. 


Figure  3.  Curves  of  changes  In  the  mean 
latitude  of  selected  European  stations. 
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Figure  4.  Curves  of  changes  of  differences  in 
the  mean  latitude  for  selected  European  stations 


Table  3.  List  of  larger  earthquakes  {M>5*5)  vicinity 


of  Belgrade. 


Nr 

Dau 

N 

E 

k 

floM) 

M 

,  ■  ,  ■  - - - - - 

Rejon  0  ) 

0  > 
1 

1964  0413 

45?3 

I8f0 

5 

! 

1 

5.5 

Skopie,  iuj. 

i  2  1 

0717 

38.0 

23.6 

ISS 

5.7 

Ort^  C'Al 

1965  0309 

39.3 

23.8 

18 

5.7 

Mone  Ejeiskie  (?J 

* 

03  31 

38.4 

22.3 

4S 

6.3 

Grec)8f‘/> 

\  s  ' 

0409 

3S.1 

24.3 

39 

5.9 

KreU  fO) 

6 

0706 

38.4 

22.4 

18 

5.8 

Grecja 

1 

1966  02  OS 

39.1 

21.7 

16 

5.6 

Crecja  C4f) 

8 

10  29 

38.9 

21.1 

1 

5.8 

Grecja  fVj 

!  9 

19670209 

39.9 

20.3 

1 

5.6 

Albania 

!  10 

0304 

39.2 

24.6 

60 

6.0 

Morse  Egejskie^^j 

1  >t 

OS  01 

39.6 

21.3 

34 

5.5 

Grecja  fV) 

i  12 

11  30 

41.4 

20.4 

21 

5.9 

Albania 

13 

19680210 

39.4 

24.9 

7 

6.0 

Morte  E*ejskie<'^^ 
Motze  Jobskiet*) 

14 

19690708 

37.5 

20.3 

0 

5.5 

!  10 

0612 

34.4 

25.0 

22 

5.8 

KxtMCkJ 

i  17 

19700408 

38.3 

226 

23 

5.8 

Orecjaf^J 

'  18 

1972  09  04 

35.2 

23.6 

14 

5.9 

KxetaC6; 

19 

0917 

38.4 

20.3 

33 

5.6 

Gnej»C¥) 

20 

0913 

38.0 

22.4 

75 

5.8 

21 

11.04 

38.9 

20.5 

13 

5.6 

Grecj»^¥) 

22 

1129 

35.2 

23.8 

37 

5.6 

Kreta  C(^,j 

Key:  (1)  Region;  (2)  Skopie,  Yugoslavia;  (3)  Yugoslavia; 

(4)  Greece;  (5)  Ae^eian  Sea;  (6)  Crete;  (7)  Albany;  (8) 
Ionian  Sea. 


Table  ^a.  Values  of  Aq  for  meridional  base 
of  Jozefoslaw  station  in  mllllgals. 


Nr 

Nt4-S14 

N6-S6 

N3-S3 

serii 

jg  1 

•HA, 

dg 

Ag 

•ha,  \ 

I 

1'“ 

-19.442  •  ±0.018 

-8.281 

±0-012 

-4.963 

1  ±0.012 

n 

-19.530  1  ±0.019 

-8.337 

±0.010 

-4.961 

±0.012 

m 

-19.441  j 

±0.029 

-8.342 

±0.016 

-4.929 

±0.006 

IV 

1  _  _ 

-19.426  j 

±0.019 

-8.497 

_ _ 

±0.016 

-4.970 

±0.010 

Table  ^b.  Values  of  q  and  Aq  for 
meridional  base  Jozefoslaw. 


Nr 

serii 

N14 

Mg 

A<p 

-41S  j  N6-S6 

tntA,  1  Mg  1  ntdA, 

mj,  j  Atf  j  /nj. 

N 

SAg 

A<p 

3-S3 

1 

iHAr  1 

ii-i  1  -"r 

1  -0;'0044 

±0.0265 

±0"0013 

-0.056 

-0:'0028 

+0.015 

±0"0007 

+0.002 

+o;'oooi 

±0.017  i 
±0:'0008  1 

—  —  [ 

!  iii-ii 

i 

+0.089 

+0;'0044 

±0.034 

±o:'ooi7 

-0.005 

-o;'ooo2 

±0.019 

±o;'ooo9 

+0.032 

+o:'ooi5 

±0.013  i 
±o;'ooo6  I 

1  iv-ni 

1 

+0.015 

+o:'ooo7 

+0.034 
j  ±0"0017 

-0.029 

-o:'ooi4 

±0.023 

±o;'ooii 

-0.041  j  ±0.012  i 

-0;'002oj  ±0"0006  j 
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Summary 

OUTereacea  of  latitudes  and  of  mean  latitudes  of  stations  located  along  a  common  meridian  were 
studied  in  order  to  find  influences  of  local  geo|)hysical  pbenomcna  on  their  variations.  Soane  conelations 
of  variations  of  the  mean  latitude  of  Belgrade  with  an  occurance  of  earthquakes  in  its  vicinity  can  be  notic¬ 
ed  u  well  as  instabilities  of  mean  latitude  of  stations  located  in  the  vicinity  of  sea  coasts. 

Results  of  eaperimental  gravimetric  meuurements  of  variations  of  a  meridian  component  of  a  vertical 
on  the  meridian  base  of  the  Jdaefoslaw  station  are  picacnted. 
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GEOMETRIC  METHOD  OF  DETERMINATION  OF  THE  SHAPE  OF  EARTH 
IN  MOUNTAINOUS  TERRAINS 

Aniela  Makowska,  Zbigniew  Zorski 
1 .  Introduction 

The  analyses  carried  out  in  works  /6/,  /?/  lead  to  the 
conclusion  that  there  is  a  possibility  of  simultaneous  determi¬ 
nation,  from  equalization  of  trigonometric  height  network,  of 
excessive  heights,  deviation  of  plumb  line  (component  deviations 
of  vertical)  and  coefficients  of  refraction  at  each  point  of 
the  network. 

This  method  of  utilization  of  the  results  of  trigonometric 
levelling  in  mountainous  terrains  provides  information  on  both 
the  building  of  atmosphere  (coefficients  of  refraction)  and 
the  weight  (mass)  field  of  Earth  (components  of  the  deviations 
of  plumb  lines)  and,  at  the  same  time,  enables  to  determine 
ellipsoidal  heights  (heights  counted  along  the  normal  to  ellipsoid 
of  the  physical  surface  of  Earth), 

In  view  of  difficulties  of  determination  of  the  shape  of 
Earth  in  mountainous  terrains  by  astronomical-geodetic  and 
gravimetric  methods, /2/,  trigonometric  levelling  may  be 
considered  as  an  auxiliary,  geometric  method  to  detennine  the 
shape  of  Earth  in  those  terrains.  However,  the  accuracy  of 
determination  of  heights  and  of  plumb  line  deviations  by  means 
of  geometric  method  is  strictly  dependent  on  accurate  knowledge 
of  the  coefficients  of  refraction. 
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The  treatment  of  the  refraction  phenomenon  in  the  form  of 
a  mathematical  model  with  constant  coefficients  of  refraction 


at  stands  (stations),  as  suggested  in  works  A/  and  /5/,  is 
a  large  step  forward  in  this  field.  However,  it  will  be  the 
subject  of  this  work  to  find  whether  the  such  accepted  mathema¬ 
tical  model  for  refraction  corresponds  to  the  real  conditions 
governing  in  mountainous  terrains.  The  analyses  will  be  carried 
out  on  the  basis  of  measurements  made  in  Polish  Tatra  Mountains, 

2,  Method  of  determination 

Measurements  of  vertical  angles  and  lengths  of  sides  were 
carried  out  in  July  and  August  1974.  The  angle  measurements  were 
done  under  direction  of  the  authoress,  and  linear  measurements 
by  Drs,  Kleczko  and  T,  Jasinski,  on  seven  points  of  the 
Tatra  network  (Figure  1),  During  the  measurements  the  weather 
was  changeable  but  mostly  sunny.  For  observation  of  angles 
we  used  theodolite  Wild  T3;  they  were  checked  before  the  use  - 
and  errors  in  the  division  of  vertical  circle  /8/  and  of 
micrometer  (of  runs  and  systematic  errors)  were  determined. 

For  precise  measurements  of  vertical  angles,  however,  special 
attention  should  be  paid  to  instrumental  and  operator  errors, 
since  in  this  type  of  instruments  there  is  no  possibility  of 
displacing  the  vertical  circle  and  micrometer  (we  measure  at 
the  same  point  of  wheel  and  micrometer,  and  we  measure  only  one 
directioi^L 


In  order  to  determine  personal  errors,  measurements  were 
made  by  several  observers  taking  part  in  the  project;  they  were 
made  simultaneously  with  several  instruments,  measuring  vertical 
angles  of  several  targets  located  at  different  heights.  Helio¬ 
tropes  and  reflectors  served  as  the  targets.  The  vertical  angles 
In  the  network  were  measured  in  four  series,  each  series  consisting 
of  four  aimings  at  two  positions  of  the  wheel  (circle).  Forty 
vertical  angles  were  measured.  The  lengths  of  sides  in  the 
network  were  measured  by  means  of  an  electrooptic  laser  range¬ 
finder  of  the  Swedish  firm  AGA  Model  8.  Standard  frequencies  of 
the  rangefinder  were  checked  before  and  after  measurements  of 
the  network.  The  value  of  the  addition  constant  and  the  linearity 
of  the  scale  of  phasometer  were  checked  on  the  standard  base  P.W. 

The  sides  were  measured  in  both  directions. 


5.  Discussion  of  results 


The  following  correction  equation  /?/  was  applied  to  each 
observed  zenithal  distance: 


•MM  sin2^fi  sm^fi  J  Sii 

Oij  2  All 


(1) 


where 


-  zenithal  distance, 

5s;  -  spatial  distance  between  points  fl.  ^ 
l^ij  -  radius  of  curvature  of  the  normal  cross-section  of 
ellipsoid  in  azimuth  •CIJ 
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In  Table  1,  Lines  1-6  contain  differences  between  ellipsoidal 

heights  calculated  on  the  basis  of  results  of  precise  levelling 

tiJt 

(orthometric  heights  ft  )  and  astronomical-gravimetric  levelling 
(interval  between  geoid  and  ellipsoid  placed  tangentially  to 
geoid  at  the  point  1-N),  and  heights  obtained  from  the  trigonomical 
levelling  of  network.  Lines  7-1 8  contain  components  of  the 
deviations  of  plumb  line  (it  was  assumed  that  at  Point  1  -fei  =  0 
and  tZi  =  0).  Lines  19-28  contain  coefficients  of  refraction 
for  each  position. 

For  Variant  I,  discrepancies  between  the  determined  heights 
are  quite  high,  up  to  13  cm.  Also,  differences  between  coeffi¬ 
cients  of  refraction  obtained  from  the  levelling  of  network 
and  from  geodetic  measurements  are  considerable,  reaching,  for 
example^  AlC  =  0,055  for  the  Point  1 ,  The  conclusion  that  follows 
is  that  the  adopted  mathematical  model  for  refraction,  which 
assumes  the  constancy  of  refraction  coefficients  at  each  position, 
corresponds  only  approximately  to  the  real  conditions  governing 
at  the  Tatra  grounds. 

We  introduced  into  equation,  therefore,  additional  unknowns 
In  view  of  differences  of  ground  cover  and  shape  for  each  line 
of  eight.  We  changed  refraction  coefficients  for  points  1  and  4 
since  the  lines  of  sight  for  these  points  were  running  partly 
over  forests  and  partly  over  sunny  Southern  slopes.  This  is  now 
Variant  II,  The  adoption  of  this  model  reduced  considerably 
discrepancies  both  for  heights  and  for  coefficients  of  refraction 
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(Table  1 ) *  The  same  nodel  was  further  Modified  by  assuaing 
only  two  constant  points,  one  with  known  height  and  deviation 
of  vertical  and  the  other  with  deviation  of  vertical 
(variant  II®,  Table  1). 

On  the  other  hand.  Variant  III  refers  to  a  nodel  with 
constant  coefficient  of  refraction  for  the  whole  network. 

Assuming  knowledge  of  height  at  Point  1  and  knowledge  of  the 
deviation  of  verticals  at  all  points,  we  determine  frMi  the 
leveling  only  differences  in  height  and  the  constant  coefficient 
of  refraction  for  the  whole  network  (Table  1), 

Summing  up  the  results  of  this  analysis  we  have  to  conclude 
that  the  acceptance  of  a  mathematical  model  for  refraction, 
with  constant  refraction  coefficients  at  a  station,  is  not 
always  appropriate,  particularly  when  the  lines  of  sight  run 
over  slopes  with  different  exposures  to  sun  and  over  terrain 
with  different  vegetation  cover  around  the  station  (points  1 
and  4),  It  is  recommended,  therefore,  to  use  in  levelling  a  model 
with  group  refraction  coefficient  at  the  station,  particularly 
for  points  lying  low  (less  than  about  1100  m  above  sea  level). 
This  procedure  provided  us  with  an  accuracy  of  determination 
of  the  components  of  deviation  of  vertical  to  ^  6®®,  and  of 
the  height  to  ^  4  cm.  The  accuracy  of  determination  of  the 
deviations  of  vertical  by  the  astronomical-gravimetric  method 
is  estimated  at  ^  3  /2/,  This  comparison  entitles  us  to  draw 
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conclusion  that  the  trigonometric  levelling  can  be  used 
as  an  auxiliary,  geometric  method  of  determination  of  the 
shape  of  Earth  in  mountainous  terrains . 


(2)hj-1013 


Figure  1. 
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Table  1.  Results  of  trigonometric  levelling  of 
elevation  network. 


Nr 

1 

II 

II* 

Ui  ■■  "1 

wariantu 

A 

"" 

m, 

A 

nil 

A 

mi  1 

dHi  (cm] 

+  8.7 

±3.3 

+  1.7 

±1.9 

+  2.2 

±3.9 

+4.0 

±3.9  i 

dHi  (cm] 

+  5.6 

±4,2 

-2.7 

±2.4 

+0,9 

±4.7 

+4.7 

±4.3 

dH*  (cm] 

-3.5 

±4,5 

-3.9 

±2,3 

-4.6 

±4,6 

0,0 

±4.2 

dHi  (cm] 

-5.6 

±3.5 

-2.6 

±1.9 

-0.4 

±2.7 

-7.5 

±3,7 

dHt  (cm] 

-9.5 

±3.5 

-2.9 

±2,0 

-2,0 

±3.6 

-11,6 

±3.9 

<W7  (cm] 

-13.3 

±3.6 

-2.7 

±2.1 

-4.1 

±3.2 

-13,9 

±4,1 

■ 

{i“> 

_ 

_ 

_ 

-2.1 

±5.3 

_ 

— 

•A**' 

~ 

- 

— 

— 

-0.2 

±3.9 

— 

— 

fV 

_ 

— 

— 

— 

+2.1 

±5.8 

— 

— 

•A"’ 

— 

_ 

-1.6 

±5.6 

— 

— 

— 

— 

— 

— 

+  7.5 

±5,3 

— 

— 

«A**‘ 

_ 

— 

— 

-2.7 

±5.7 

— 

— 

{I..1 

— 

— 

~ 

+24,1 

±3.4 

— 

— 

•A“> 

— 

— 

+  10,5 

±3.4 

— 

— 

{i"> 

-- 

— 

— 

— 

— 

— 

— 

— 

uk 

— 

— 

_ 

+  10,5 

±5.5 

— 

«A"> 

— 

— 

+  1.1 

±4.2 

- 

1 

JiC. 

1 

1 

0.159 

±0.006 

0,155 

±0,010 

JIT', 

0.125 

±0.006 

0,1% 

±0,004 

0,120 

±0,006 

AK\' 

I 

1 

0,096 

±0,006 

0,093 

±0,008 

AKi 

0.125 

±0,005 

0,118 

±0,003 

0,112 

±0,010 

AKi 

0,117 

±0,006 

0.107 

±0,004 

0,107 

±0,011 

AK» 

) 

0,097 

±0,003 

0,107 

±0.011 

AKi 

1  0,091 

J  ±0,005 

0,114 

±0,004 

0,127 

±0.013 

AKs 

0,127 

±0.006 

0,133 

±0.004 

0,137 

±0,006 

AK, 

0.126 

±0,005 

0,134 

±0,003 

0,133 

±0,005 

AKi 

0,120 

±0,004 

0,126 

±0,002 

0,120 

±0,009 

;l 

1 

me 

±5 

•'.2 

±2”  .7 

±2“  ,6 

±7 

•< 

.9 

29 
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Summary  | 

In  the  article  the  refraction  formula  for  trigonometric  levelling  in  mountains  is  analysed. 

It  it  shown  in  the  work  that  constant  coefficient  of  refraction  is  not  always  appriopriate  in  survey,  ' 
particularly  when  lines  of  sight  run  over  the  sun  heated  slopes,  or  over  the  ground  around  a  statioD  with  j 
different  vegetation  cover.  In  this  case  the  use  of  group  refraction  coefficients  is  suggested  in  particular 
for  points  situated  below  1100  m  above  sea  level.  This  analysis  is  based  on  observations  carried  out  on 
the  Pnliik  region  of  Tatra  MounUins. 


-30- 


INFLUENCE  OF  THE  REFRACTION  COEFFICIENT  OF  ELECTROMAGMETIC  WAVES 


ON  RESULTS  OF  MEASUREMENTS  WITH  ELECTROMAGNETIC  RANGEFINDERS 

Michal  Knoll 

1 .  Introduction 

The  phenomenon  of  refraction  is  important  in  geodesy 
whenever  there  is  the  case  of  precise  measurements  carried 
out  with  optical  instruments  (theodolit,  leveller,  electro¬ 
optic  rangefinder). 

Many  authors  dealt  with  the  problem  of  refraction  considering 
it  from  the  viewpoint  of  geometrical  optics.  Such  an  approach 
is  right  if  it  can  be  assumed  that  the  wavelength  =  0, 

This  assumption  can  be  accepted  in  the  range  of  visual 
radiation  X  =  0,4-0,76  ^m.  But  we  cannot  make  such  an 
approximation  for  microwave  rangefinders  workinv  on  centimeter 
waves , 

When  a  side  is  measured  with  two  instruments  we  shall  get, 
on  the  whole,  two  results.  The  ratio  of  these  values  is  called 
the  scale  of  the  side. 

In  the  last  years,  many  measurements  have  been  made  using 
the  microwave  and  optical  rangefinders.  The  comparison  of 
i^esults  between  them  reveals  the  occurrence  of  a  systematic 
error  of  the  scale.  In  the  majority  of  cases,  measurements  in 
the  optical  spectrum  give  higher  values  of  distance  than 
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meaeurements  in  microwave  range.  The  work  /8/  gives  results  of 
measurements  of  two  the  same  sides  on  two  Polish  experimental 
networks  performed  with  the  microwave  and  electrooptical 
rangefinders,  A  hypothesis  is  put  forth  about  existence  of 
a  difference  of  scales.  As  a  result  of  the  analysis  carried  out 
it  was  established  with  98^  probability  that  there  are  no  reasons 
to  reject  the  accepted  hypothesis. 

An  incorrect  determination  of  distance  with  a  rangefinder 
may  be  caused  by;  instability  of  electronic  system,  careless 
measurement,  insufficient  knowledge  of  the  coefficient  of 
refraction,  and  inappropriate  comparation,  A  unit  of  length 
vhich  is  conferred  upon  the  rangefinder  during  the  comparation 
is  transferred  onto  sides  of  triangulation  network.  It  is  assumed 
that  this  unit  is  preserved  by  electronic  system  of  the  instrument 
rlth  accuracy  not  lower  than  10*^,  If  we  assume,  moreover,  that 
the  electronic  system  of  rangefinder  is  efficient,  and  that  the 
comparation  and  measurement  were  carried  out  properly,  the  only 
source  of  errors  causing  a  systematic  scale  difference  can  be 
the  refraction  coefficient  of  electromagnetic  waves  determined 
with  Insufficient  accuracy.  We  shall  consider,  therefore, 
the  effect  of  various  factors  on  values  of  refraction  coeffi¬ 
cient  for  light  waves  and  microwaves. 
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?.  Refraction  coefficient  In  a  dleperelon  medium  for  light  wares 

2,1 ,  Refraction  coefficient  for  standard  air 

A  formula  describing  the  coefficient  of  refraction  for 
standard  air  was  introduced  in  1 953  by  B,  Edlen  on  the  basis  of 
the  equation  of  Cauchy  and  studies  of  Barrell  and  Sears  A/. 
Subsequently,  studies  in  this  field  were  published  by  Rank  (1958), 
Svensson  (I960),  Peck  (1962)  and  Erickson  (1961),  providing  new 
information  about  dispersion  of  electromagnetic  waves  in  air 
and  in  water  vapor. 

On  the  basis  of  the  latest  studies  /5/  the  dispersion  formula 
for  standard  air  has  the  form: 


(n-l)lO*  =  8342,13  +  2406030(130-(T*)‘'  +  l5997(38,9-<r’)- 


-1 


(2.1) 


where  ^  -  wavenumber  in  vacuo  in  ;im"’  , 

2,2.  Dependence  of  the  refraction  coefficient  Iv 
on  temperature  t  and  pressure  ^ 

B,  Edlen  in  reference  /5/  introduces  formula  which  gives 
dependence  of  Vv  on  and  p  .He  starts  from  the  relation 


(*•— IXr  “  KxOfp, 


(2.2) 


where  -  a  factor  dependent  on  wavelength,  and  -  a  density 
factor' dependent  on  temperature  and  pressure, 

f . 

The : proof  is  based  on  the  basic  Lorents-Lorenz  relation.  The 
equation  has  the  form; 
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l+j»(0.817-0,0133/)10-‘ 


720,775 


l+0,0036610r 


(2.5) 


and  can  be  substituted  with  the  fonaula  of  Barrell  and  Sears 

o 

of  1939*  with  accuracy  of  the  determination  of  n  of  1  x  10  : 


(n— i)„  =  (n— 1), 


0,001 38823/7 
1+0,003671/ 


(2.4) 


The  substitution  of  formula  (2.5)  with  fenewla  (2.4)  with  the 
above  given  accuracy  is  valid  if  lies  within  the  range 

5-30  °C,  and  ^  in  the  range  700-800  Tr, 


2,3.  Effect  of  the  content  of  water  vapor 


Following  Barrell  and  Sears  /!/  we  shall  write: 


.  Kx{p-f)[X+P,(p-f)]  .  +?',/) 

'  1  +  a/  1  +  */ 


(2.5) 


This  equation  can  be  transformed  /8/,  by  neglecting  expressions 
as  values  of  the  order  of  10”^  at  the  most,  in  order 
to  obtain  members  dependent  on  humidity  ^  and  members  inde¬ 
pendent  on  ^  ,  We  shall  get  then  an  equation  which  provides 

effect  of  taking  account  of  humidity  on  the  accuracy  of  determi¬ 


nation  of  refraction  coefficient: 


(2.6) 


Values  and  ^  refer  to  dry  air,  and  values  fC*^  and  - 
to  moist  air.  They  are  defined  as  follows: 
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ATt 


(•-U 


0  +  15.), 


760 

(31,59  +  0, 2963(t*)10-*. 
(0,817-0,0133/)  10-‘, 
27xl0-‘, 

1 


273,15 


Equation  (2,6)  is  well  satisfied  in  the  limits  of  spectrum 
644.0-435.9  nm. 

Recent  studies  of  Erickson  /6/  determined  refraction  of  water 
vapor  with  ^rreater  accuracy  and  for  a  larger  range  of  spectrum 
644.0-361 .1  nm.  The  equation  as  a  function  of  wavelength  has 
the  form; 


n„f-n„  =  -/(5,722  -0,0457<T»)10-*  (2,7) 

and  is  valid  for  conditions  which  do  not  diverge  much  from 
t  =  20  °C,  p  =  760  Tr,  1=10  Tr. 

2,4.  Effect  of  the  content  of  carbon  dioxide 

It  is  assumed  that  in  the  standard  air  the  content  of  CO2 
is  0,0003  parts  per  unit  of  volume.  This  corresponds  to  partial 
pressure  0,23  Tr, 

B,  Edlen  /4/  gives  formula  by  means  of  which  he  calculated 
the  refraction  coefficient  of  air,  containing  x  parts  CO^ 
per  unit  volume  in  relation  to  the  standard  air; 

(.-1),  .  (1  +0,540(jf-0,0003))(«-l)b.  (2,8) 
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The  difference  between  refraction  coefficients  for  CO^  and 
for  the  standard  air,  as  a  function  of  wavelength,  has  the  fora: 

(14485 +  n7,8o»)10-*  (2,9) 

Values  calculated  from  the  above  formula,  given  by  Edlen,  agree 
with  values  reported  by  T,  Masui  /9/, 

3.  Refraction  coefficient  of  microwaves 

The  vector  of  electrical  displacement  is  connected  with 
the  vector  of  potential  ^  in  a  given  medium  through  the  formula; 

/)  =  «.£.  (3.0 

where  -  dielectric  permeability  of  the  medium. 

In  isotropic  medium  we  can  write; 

a-e.£.  (5.2) 

Making  use  of  the  dielectric  permeability  of  vacuum  Sa 
permeability  of  air  ,  the  equation  (3.2)  can  be  written; 


Let  us  denote; 


D  >  t,E  »  toE+tfE, 


*0 
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The  relative  permeability  of  air  can  be  written  as  the  sum  of 
permeabilities  of  dry  air  and  water  vapor  : 


<- (3*3) 

where  A/  “  number  of  particles  in  unit  of  volume. 

Particles  of  water  vapor  have  a  constant  dipole  moment. 

Its  effect  on  the  field  will  depend  on  temperature  .  Moreover, 
an  induction  moment  appears  in  the  electromagnetic  field.  If  we 
take  this  into  account  in  (3«3)  end  if  we  substitute  appropriate 
values  for  Nj  and  Nw  ,  we  get; 


where  a  constaint  connected  with  induction  moment,  - 

—  a  constant  connected  with  dipole  moment,  p  -  pressure  of 
dry  air,  -  partial  pressure  of  water  vapor.  After  transfor¬ 
mation  of  (3.4)  we  obtain: 


»  =  I  +  e.*- +  4 -4- 

T  T  T* 

Utilizing  the  relation: 


(3.5) 


we  obtain; 


.  e— 1 
1  +  — —  + ... 


^  I  e,a-‘  p  e„a-‘  /  f 

2  2  r  2  r*  • 


(3.6) 
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For  radio  waves  we  use  usually  a  value  called  refraction 
indicator: 

N-  (n-DlO*. 

Introducing  into  (3.6)  the  constants: 

j— 

we  obtain: 

N’‘A^+b^+c-L.  (3,7) 

In  practice  we  meet  often  the  relation: 

(3,8) 

Constants  A .  B  ,  C  are  determined  experimentally. 

Formulas  (3.7)  and  (3.8)  have  not  been  deduced  strictly. 

Most  of  all,  the  calculations  do  not  take  into  account  the  fact 
that  the  dielectric  permeability  is  a  function  of  the  frequency 
of  field. 

k.  Analysis  of  equations  describing  the  refraction  coefficient 
4.1 .  Formula  for  microwaves 

Differentiating  equation  (3.7)  with  respect  to  T*  ,  |^  and  "f 
and  substituting  the  values  A  =  103.49»  h  =  -17.23,  C  =  495jB22, 
for  conditions  J  =  25  Tr,  ^  =  760  Tr,  "7^  =  303  we  obtain: 


-■*  i 
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-AT-*p-BT-^f+CT-*f=  -1,7432, 


£L 

dT 

dN 

dp 


*  AT-' 


0,3415, 


dN_ 


BT-'  +  CT-*  *  3,3437. 


Assuming  errors  in  determination  of  the  parameters  of  atmosphere  to  be 
/T  =  1  °K,  oCp  =  1  Tr,  =  1  Tr,  we  obtain:  I 

=  5.631  X  10“^, 

which  gives  the  relative  accuracy  of  determination  of  the 
length  of  side  1:l7';j000. 

4.2.  Formula  for  light  waves 

Differentiating  equation  (2,5)  with  respect  to  t"  $  f  and  ^ 
and  substituting  suitable  values  for  conditions  ^  =  25  Tr, 

^  =  760  Tr,  ^  =  30  °C  we  obtain: 


dn 

dp 

dn 

If 

dn 

IT 


(1  +2/;A -.2/A)  «  3,49 x}0‘\ 

fC*  IT 

(1+2//?!)+  (2/(J.-l  -2pW  =  -7.54X  10-*, 

1  +  af  l+a^ 

{x,»(;>-/)li+A(/>-/)]+«/(i+^y^}~^—  =  8,69x  10-’. 

(1+ar)* 
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Aseunlng  error?  in  determination  of  the  parameters  of  atmosphere 
to  be  J(  °c,  s  1  Ti* f  ^  Tr f  W0  ob^Aliie 

4>V  =  0.94  X  10"^ 

which  gives  the  relative  accuracy  of  determination  of  the 
length  of  side  1;1,064»000, 

Table  1  compares  relative  errors  of  the  measurement  of  the 
side  with  microwave  and  light f (optical)  rangefinders  depending 
on  values  dp  .  4  and  4i  . 

In  order  to  determine  error  caused  by  neglecting  humidity 
irhen  calculating  the  refraction  coefficient  for  light  waves 
fre  shall  make  use  of  the  formula  (2,6),  For  “t  =  30  p*  =  25  Tr, 
^  =  760  Tr  we  obtain: 

Hiff  -  =  -  1,54  X  10"^ 

^om  Erickson's  formula  (2,7)  we  have: 

The  effect  of  CO^  on  the  refraction  coefficient  can  be 
determined  on  the  basis  of  equations  (2,8)  or  (2,9).  If  we 
assume  ten  times  larger  concentration  of  CO^  than  in  the 
standard  air,  which  is  a  very  large  value,  both  equations  give 
numbers  significant  only  at  the  seventh  decimal  point. 
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Concluding  reaarke 


The  dispersion  of  electromagnetic  waves  In  atmosphere 

af 

Is  caused  by  absorption  and  polarization  ^ ^molecules. 

Neither  polarization  nor  dielectric  constant  depend  on 

Q 

frequency  if  i  It  does  not  exceed  the  value  of  10  Hz. 

A  sharp  drop  of  value  is  observed  in  microwave  region 
(  U  =  10^  -  10^^  Hz). 

For  the  wave  frequency  bd  =  10^^  Hz  and  higher  -  hence 
for  light  waves  -  the  orientational  polarization  (connected 
with  particles  possessing  permanent  dipole  moments)  disappears, 
and  also  induced  polarization  partly  disappears,  and  it  is  only 
electron  polarization  that  remains /10/.  The  value  of 
does  not  differ  much  from  that  of  molar  refraction,  which  can  be 
calculated  from  the  Lorentz-Lorenz  equation.  For  microwaves, 

therefore,  this  part  of  refraction  will  be  important  which  is 
connected  with  orientational  polarizability. 

Van  Vleck  /I  2/, /I  3/  foresaw  two  absorption  spectra  for 
microwaves.  One  can  observe  these  spectra  for  wavelength  1 .35  cm 
/II/  (lines  of  water  vapor  between  two  rotational  states)  and 
wavelengths  4  mm  and  6  mm  /2/  (absorption  of  magnetic  dipoles 
of  0^  molecules). 

Lately,  carrier  frequencies  of  microwave  rangefinders 
increased  considerably  and  approached  the  resonance  frequency 
of  water  vapor  (  X  =  1.35  cm);  it  can  be  assumed  then  that 
disturbances  in  dispersion  are  caused  by  phenomenon  of 
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anonalous  dispersion*  This,  in  turn,  may  cause  changes  in 
frequency  of  the  emitted  and  received  waves  by  the  rangefinder 


(Figure  1).  A  change  of  frequency  means  a  shift  of  the  maximum 
energy  ^  into  another  part  of  the  spectrum. 

It  appears  desirable,  therefore,  to  determine  dlspersional 
formulas  for  each  model  of  the  rangefinder, 

A  large  influence  upon  the  measurement  results  with  microwave 
rangefinder  can  be  exerted  by  the  so  called  near-ground  effect. 
This  effect  is  a  nonlinear  dependence  of  the  pressure  of  water 
vapor  on  height  near  the  surface  of  earth.  On  the  whole,  the 
line  of  sight  runs  higher  over  the  terrain  in  its  middle  part 
than  at  the  end  points.  The  averaging  of  humidity  measured  at 
the  end  points  does  not  appear,  therefore,  justified  /?/. 

Measurements  in  the  microwave  region  in  industrial  centers 
should  be  avoided.  The  pollution  of  atmosphere  with  chemical 
compounds,  occurring  there,  may  lead  to  anomalous  dispersion. 

It  would  appear  also  that  the  results  of  microwave 
measurements  could  be  improved  if  the  relation  t*  cf«) 
was  taken  into  account  in  formula  (3.6).  This  relation  could, 
in, particular,  change  the  value  in  equation  (3.4). 
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Table  1 


tTr) 

(TrI 

a 

[deg] 

Didmien  twktluy 

dt 

s 

DabnicR  mikrafalowy 

ds 

s 

0.3 

0,3 

0,3 

1:3  364  000 

1:561  000 

1 

1 

1 

1:1  064  000 

1:177  000 

1  1 

2 

2 

1:561000 

. 

1:89  000 

• 

4 

4 

1:285  000 

1:44  000 

2 

5 

5 

1:143  000 

]:36000 

_  .  1 

Key:  (1)  Light  (optical)  rangefinder;  (2)  Microwave 

rangefinder. 

i  »  _  ’ 


I  /' 


\  '  ! 
\  1  \2 
\  I ' 

Vjiy 


frequency  sent  by  rangefinder 
(max,  E  for  ) 
resonance  frequency 
""  (min.  E  for  ttr  ) 

frequency  received 

(max.  E  for  ) 


Figure  1 
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Summary 

There  is  an  attempt  of  explanation  of  a  phenomenon  that  exists  in  the  form  of  systematic  diflerence 
of  scales  consequent  on  microwave  and  electrooptical  distance  measurements.  There  is  assumed  that  the 
essential  factor  causing  this  phenomenon  is  an  allowance  for  refraction  coefficient  with  insufficient  accuracy. 
There  is  discussed  the  formulae  for  refraction  coefficient  in  atmosphere  in  the  two  spectrum  regions.  For 
light  waves  there  are  considered  the  formulae  by  Edlen  and  Erickson  and  the  formula  by  Bareli  and  Sears 
that  is  so  far  applied  to  geodesy.  There  are  pointed  out  some  inaccuracies  being  consequent  on  deriving 
of  the  formulae  for  microwaves.  There  is  also  discussed  the  phenomenon  of  absorption  in  microwave 
region.  It  is  proved  that  absorption  has  the  essential  influence  on  the  systematic  difference  of  scales. 
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INVESTIGATION  OF  THE  INFLUENCE  OF  ERRORS  OF  JUNCTION  POIMTS 


£ 


r>W  THE  ACCTOACY  OF  ADJDSTED  GEODETIC  CONSTRUCTIONS 
Jerzy  Mlalduo,  Kazlmlerz  Slkorskl 

1 .  Introduction 

When  adjusting  geodetic  networks  It  Is  assumed  as  a  rule 
that  the  junction  points  are  exact,  without  errors.  This  Is 
dictated  by  the  necessity  of  maintaining  Invariable  coordinates 
In  the  existing  catalogs.  The  necessity  arises  mainly  from  the 
adopted  system  of  collecting  and  storing  geodetic  data. 

At  present,  as  the  measurement  techniques  and  the  methods 
of  treatment  of  measurement  results  continue  to  be  developed. 

It  occurs  more  and  more  frequently  that  networks  of  the  lower 
orders  are  'characterized  by  a  higher  Internal  accord  than  the 
supportive  networks.  The  attachment  of  these  lower-order  sets 
to  networks  characterized  by  lower  accuracy,  with  simultaneous 
assumption  of  the  errorless  nature  of  the  points  in  connected 
networks,  results  in  a  considerable  deformation  of  both  the 
adjusted  measurement  results  and  coordinates  of  the  new  points, 
and  of  the  estimate  of  their  accuracy.  The  problem  of  adjustment 
with  consideration  of  the  effect  of  errors  of  junction  points 
was  treated  among  others  In  the  works  /I/  ap^^4/. 

An  Interesting  problem  was  considered  in  work  /2/,  The 
author  studied  the  effect  of  the  errors  of  Junction  points 


and  of  their  placenent  on  the  accuracy  of  the  nevly  detemined 
points.  He  coses  to  the  conclusion  that  when  adjusting  the 
nunher  of  points  in  newly  established  network  it  is  necessary 
to  limit  the  nuJBber  to  absolute  necessity  (about  3)* 

In  the  present  work,  studies  on  the  effect  of  the  nuaber 
of  junction  points  of  connections  are  carried  out  under  the 
assumption  of  the  presence  of  errors  and  of  the  dependence 
between  these  points. 

2.  Construction  of  supportive  network 

In  order  to  investigate  the  effect  of  connective  points 
on  errors  of  new  established  points  we  constructed  a  supportive 
network  shown  in  Figure  1 ,  The  same  Figure  contains  also  ellipses 
of  errors  of  the  points  in  supportive  network.  In  the  construction 
of  network  we  were  guided  by  the  following  premises: 

(a)  the  determined  points  are  placed  symmetrically  with 
respect  to  poihts  A,  B  which  define  the  system, 

(b)  points  of  the  supportive  network  are  also  symmetrical 
with  respect  to  the  points  of  connected  network  (consisting  of 
four  points  Nos.  1,  2,  3»  k), 

(c)  all  the  constructions,  both  in  the  supportive  and  in 
connected  network!  consimt  of  isosceles  triangles, 

(d)  rejection  of  each  pedr  of  e3rmi^dtrlcally  placed  junction 
points  causes  reduction  of  the  nuaber  of  measurement  results 

by  the  same  value. 


The  sane  accuracy  of  angular  measurements  was  adopted  In 
adjustment  of  supportive  network  and  connected  network.  At  the 
same  time  it  was  assumed  that  angular  measurements  were  done 
from  each  determined  (defined)  point  to  all  connective  points 
lying  on  the  same  side  of  the  axis  of  symmetry  as  the  determined 
(defined)  point, 

3.  Adjustment  of  supportive  network  and  connected  network 

The  supportive  network  was  adjusted  by  an  indirect  method, 
determining  the  matrix  of  covariance  of  point  coordinates 
in  the  form 

Cav(jr).  (A' A)-'. 

It  was  assumed  during  the  adjusting  that  the  elements  of  both 
supportive  and  connected  network  were  determined  with  the  same 
accuracy. 

Investigation  of  the  effedt  of  errors  was  carried  out  in 
two  groups  A  and  B  and  six  variants  in  each  group.  The  supportive 
(determined)  points  in  each  group  were  the  same. 

In  the  group  of  variants  A  it  was  assumed  that  the  junction 
points  were  errorless  and  independent  of  each  other;  in  the 
group  B  -  that  there  were  errors  and  mutual  dependence. 

Two  or  four  pairs  of  points  lying  symmetrically  with  respect 
to  the  connected  network  were  rejected  (omitted)  in  each  variant, 
causing  the  same  change  in  error  equations  in  appropriate  variants. 
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For  junction  points  in  both  groups  of  variants  we  took: 


yariant  1  101, 102, 103, 104 

Variant  2  101, 102, 103, 104, 105, 106, 107, 108 

Variant  3  101, 102, 103, 104, 105, 106, 107, 108, 109,  110,  111,  112 

Variant  4  105,  106, 107, 108 

Variant  5  105,  106, 107,  108, 109,  110,  111,  1 12 

Variant  6  109,110,111,112 


The  adjustment  and  the  analysis  of  accuracy  were  perforaed 
on  the  basis  of  an  algorithm  presented  in  /3/»  determining 
values  of  the  weight  coefficients  of  the  points  of  network. 
Moreover,  the  following  values  were  also  determined; 


sp(0-‘)  for 


«r- 


/Sp(g-‘) 

H 


where  -  average  value  of  the  weight  coefficient  of 

a  mean  arbitrary  point  of  network,  tX?  ~  covariance  matrix, 
-  number  of  points  of  the  network. 

The  results  of  adjustment  of  the  supportive  network  are 
presented  in  Table  1 , 
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Conclusions 


The  analysis  of  adjustment  results  allows  us  to  state  that 
in  the  considered  network,  both  in  the  variant  group  A  and  B, 
an  increase  of  the  number  of  supportive  points  improves  only 
little  the  accuracy  of  placement  of  the  new  points. 

This  conclusion  agrees  with  that  in  the  work  /2/, 

The  attachment  of  a  new  network  to  supportive  points, 
assumed  to  be  errorless  (variant  A) ,  lying  at  a  considerable 
distance  from  the  new  points  and,  at  the  same  time,  characterized 
by  a  low  accuracy,  worsens  considerably  the  accuracy  of  placement 
of  the  new  connected  points. 

In  the  case  when  we  do  not  assume  the  errorless  character 
of  the  connected  points  (variant  B) ,  the  accuracy  of  new  points 
also  worsens  but  to  a  considerably  smaller  degree. 

As  a  result  of  joint  adjustment  of  coordinates  of  the  points  of 
supportive  and  connected  networks,  one  obtains  an  improvement 
in  the  accuracy  of  location  ■^f  the  points  in  supportive  network 
by  about  10^,  on  the  average. 

The  considered  case,  however,  deals  with  too  small  a  network 
to  justify  generalization  of  conclusions.  Further  investigations 
are  required,  and  they  will  be  carried  out  by  the  authors. 
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Table  1.  Supportive  network 


Nr 

punktu 

V9„ 

101 

0.145 

0,145 

0,205 

102 

0,103 

0,103 

0,145 

103 

0,145 

0,145 

0,205 

104 

0,103 

0,103 

0,145 

105 

o,m 

0,308 

0,435 

106 

0,229 

0,229 

0,324 

107 

0,308 

0,308 

0,435 

108 

0,229 

0,229 

0,324 

109 

0,513 

0,513 

0,726 

no 

0,410 

0,410 

0,580 

111 

0,513 

0,513 

0,726 

112 

0,410 

0,410 

0,580 

S|>(Q-‘)  -  2,441,  -  0,451. 


Results  of  adjustment  of  the  connected  network 
are  shown  In  Tables  2,  3> 


Table  2.  Adjustment  results  of  a  new  network,  with  the 
assumption  of  errorles  character  and  independence  of 
coordinates  of  the  points  in  supportive  network. _ 


Nr 

punktu 

WARIANT 

A1 

WARIANT 

A2 

WARIANT 

A3 

WARIANT 

A4 

WARIANT 

A5 

WARIANT 

A6 

|/9»x+9»» 

V  9«+9w 

1 

0,145 

0,139 

0,139 

0,833 

0,790 

2,766 

2 

0,188 

0,179 

0,177 

0,927 

0,872 

2,913 

3 

0,145 

0,139 

0,139 

0,833 

0,790 

2,766 

4 

0,188 

0,179 

0,177 

0,924 

0,872 

2,913 

Table  3.  Adjustment  results  of  a  new  network,  with  the 
assumption  of  the  presence  of  errors  and  of  dependence 
of  coordinates  of  the  points  In  supportive  network. 


Nr 

WARIANT 

B1 

WARIANT 

B2 

WARIANT 

B3 

WARIANT 

B4 

WARUNT 

B3 

WARIANT 

B6 

punktu 

}/g..+Qn 

^lMM+4n 

V«MM+gn 

ifimm+Qn 

1 

0173 

0.169 

0.168 

0.236 

0,233 

0.372 

2 

0.123 

0.207 

0.206 

0.280 

0.274 

0.418 

3 

0.173 

0.169 

0.168 

0,236 

0.233 

0.372 

4 

0.213 

0.207 

0.206 

0.280 

0.274 

0.418 

101 

0.201 

0.201 

0.200 

102 

0.141 

0.141 

0.141 

103 

0.201 

0.201 

0.200 

104 

0.141 

0.I4I 

0.141 

103 

0.417 

0.416 

0,417 

0,413 

106 

0.314 

0.314 

0.306 

0,303 

107 

0.417 

0,416 

0.417 

0.413 

108 

0.314 

0,314 

0.306 

0,303 

109 

0,680 

0.699 

0,643 

no 

0,330 

0.338 

0,308 

111 

0,680 

0.699 

0,643 

112 

0,330 

0.338 

0,308 

Table  4.  Average  value  of  the  mean  arbitrary  weight 
coefficient  of  the  network  q  and  trace  of  matrix  of 
the  adjusted  system. 


WARIANT 

1 

2 

3 

4 

3 

6 

cn 

SM  dowivorwaiiia  «i, 
po  wyrOwnaniu  Sp(Q~') 

A 

0.>68 
0.1 13 

0,160 

0,103 

0,139 

0,101 

0.881 

3,106 

0,832 

2,769 

2,840 

32,273 

SM  dowi^aqrwuia  fb 
po  wyiOwmaniu  Sp(0~') 

B 

0.I9S 

(M32 

0,189 

0,143 

0,188 

0,141 

0,239 

R268 

0,234 

0039 

0.323 

0,423 

Sie6  atodclowa 

wyiOwnanicoi  Sp(fi~') 

- 

A 

0,178 

0,126 

0,299 

0^714 

0,431 

2.441 

0^383 

0,388 

0,338 

2,313 

0,637 

1,727 

Sie6  moddowa  «b 

po  wyrOwnaiiiu  SpC2~') 

B 

0,174 

0,121 

0,288 

0,666 

0,427 

2,193 

0,366 

0,333 

0,316 

2,130 

0,379 

1,343 

Key:  (1)  Connected  network  after  adjustment;  (2) 

Model  network  before  adjustment;  (3)  Model  network 
after  adjustment. 
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Figure  1 
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Summary 

Thlt  paper  djacusMa  an  influence  of  junction  pointa  errors  on  the  accuracy  of  connected  network. 
Inveatigationa  have  been  carried  out  in  the  two  groups  of  variants;  with  the  assumption  of  errorlessness 
of  junction  pointa  (A  variants)  and  with  allowance  for  erroneousness  and  dependence  of  those  points 
(B  variants). 

Resulting  Iheae  analytes  it  has  been  determined,  that  increasing  of  a  number  of  junction  points  in- 
fluences  very  small  on  increasing  of  the  accuracy  of  adjusted  network  and  that  junction  to  distant  points 
makes  worse  the  accuracy  of  connected  network  as  well. 
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